Limits on spin-dependent WIMP-nucleon cross-sections from 
the XENONIO experiment 
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XENONIO is an experiment to directly detect weakly interacting massive particle (WIMPs), which 
may comprise the bulk of the non-baryonic dark matter in our Universe. We report new results for 
spin-dependent WIMP-nucleon interactions with ^^^Xe and ^^^Xe from 58.6 live-days of operation 
at the Laboratori Nazionali del Gran Sasso (LNGS). Based on the non-observation of a WIMP signal 
in 5.4 kg of fiducial liquid xenon mass, we exclude previously unexplored regions in the theoretically 
allowed parameter space for neutralinos. We also exclude a heavy Majorana neutrino with a mass 
in the range of ^10 GeV/c'^-2 TeV/c^ as a dark matter candidate under standard assumptions for 
its density and distribution in the galactic halo. 
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PACS numbers: 95.35.-|-d, 29.40.Mc, 95.55.Vj 

Evidence for a significant cold dark matter component 
in our universe is stronger than ever [J, [1, a well- 
motivated particle candidate being the lightest neutralino 
from super-symmetric extensions to the Standard Model 
[3| • Such a particle is neutral, non-relativistic, stable, and 
more generally classified as a Weakly Interacting Mas- 
sive Particle (WIMP). The open question of the nature 
of WIMPs is being addressed by numerous direct and 
indirect detection experiments [a d, 0] ■ 

Among these, the XENONIO experiment aims to 
directly detect galactic WIMPs scattering elastically 
from Xe atoms. Moving with velocities around 10~^c, 
WIMPs can couple to nucleons via both spin-independent 
and spin-dependent (axial vector) interactions. Spin- 
independent WIMP-nucleon couplings are in general 
smaller than axial vector couplings 4]. However, for 
low momentum transfer, they benefit from coherence 
across the nucleus, and therefore the overall event rate for 
WIMP interactions is expected to be dominated by the 
spin-independent coupling for target nuclei with A>30. 
The sensitivity of XENONIO to spin-independent inter- 
actions is published in Q. 

We report here on a spin-dependent analysis of 58.6 
live-days of WIMP-search data, taken in low-background 
conditions at LNGS, which provides ~3100 meters water 
equivalent rock overburden. XENONIO is a dual phase 



(liquid and gas) xenon time projection chamber, discrim- 
inating between the predominantly electron-recoil back- 
ground and the expected nuclear-recoil WIMP signal via 
the distinct ratio of ionization to scintillation for each 
type of interaction Q . A nuclear recoil energy threshold 
of 4.5 keV was achieved, and 10 candidate events were 
recorded for an exposure of about 136 kg days after anal- 
ysis cuts (the fiducial mass was 5.4 kg). Although all ob- 
served events are consistent with expected background 
from electron recoils (see [1| for details on the analysis 
and the candidate events) , no background subtraction is 
employed for calculating the WIMP upper limits. In the 
following analysis, we use identical data quality, fiducial 
volume, and physics cuts as reported in [8|. 

For axial WIMP-nuclei interactions, the WIMPs cou- 
ple to the spins of the nucleons. Although the inter- 
action with the nucleus is coherent (as it is in the spin- 
independent case) in the sense that scattering amplitudes 
are summed over nucleons, the strength of the interaction 
vanishes for paired nucleons in the same energy state. 
Thus only nuclei with an odd number of nucleons will 
yield a significant sensitivity to axial WIMP-nuclei inter- 
actions [9] . Almost half of naturally occurring xenon has 
non-zero nuclear spin: ^^^Xe (spin-1/2) makes up 26.4%, 
and i3ixe (spin-3/2) another 21.2%. The differential 
WIMP-nucleus cross section for the spin-dependent in- 
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teraction can be written as 



da 



C. 
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where Gp is the Fermi constant, v is the WIMP velocity 
relative to the target, 5'(|q|) is the spin structure function 
for momentum transfer q > 0, and S(0) represents the 
zero-momentum transfer limit. The so called enhance- 
ment factor, C'spin, is given by 



[ap(Sp) + an(Sn)] 



,J + 1 
J ' 



(2) 



where J is the total nuclear spin, ap and a„ are the ef- 
fective WIMP-nuclcon couplings (these depend on the 
quark spin distribution within the nucleons and on the 
WIMP type), and {Sp,n) = {N\Sp.n\N) are the expecta- 
tion values of the spin content of the proton and neutron 
groups within the nucleus. 

Detailed nuclear shell-model calculations for (Sp^n) for 
two different Hamiltonians describing the nuclei exist in 
the literature ■ The Hamiltonians are based on realis- 
tic nucleon-nucleon potentials, Bonn A [ll[ and Nijmegen 
II [l^. The accuracy of these calculations is typically 
assessed by comparing their predictions of the magnetic 
moment fi with experimental values, due to the similarity 
between the magnetic moment operator and the matrix 
element for WIMP-nucleon coupling. 

For ^^^Xe the magnetic moment agrees to within 25% 
and 11% for the Bonn A and Nijmegen II potentials, re- 
spectively, using the standard free-particle g-factor (and 
to within 19% and 52% for using effective g-factors). We 
calculate our main WIMP-nucleon exclusion limits us- 
ing the Bonn A potential (Fig. [U solid curves), with 
(Sp) = 0.028 and (S„) = 0.359. In order to indicate 
the level of systematic uncertainty associated with the 
different models we also calculate limits using the al- 
ternate Nijmegen II potential (Fig. [1] dashed curves), 
with (Sp) = 0.0128 and (S„) = 0.300. For ^^iXe, the 
same Bonn A and Nijmegen II models predict the mag- 
netic moment to within 8% and 50% of the measured 
value, respectively. However, for this isotope there are 
calculations by Engel using th e Q uasiparticle Tamm- 
Dancoff approximation (QTDA) jl3| , which yield a mag- 
netic moment within 1% of the experimental value. We 
follow [iS[il and use the calculation by Engel, choosing 



-0.041 and (S^ 



-0.236, the effect of the 3 



different models for ^"^^Xe on the variation in the exclu- 
sion limits being quite small. We also note that recent 
calculations by Kortelainen, Toivanen and Toivanen [l^ 
yield values for the magnetic moments within about 20% 
and 10% for ^^^Xe and ^'^^Xe, respectively, without using 
effective g-factors. For the sake of brevity and compar- 
ison with other experimental results (which follow 
we use the (Sp^n) values detailed above. 

In the limit of zero momentum transfer the WIMP 
essentially interacts with the entire nucleus. Once the 



momentum transfer q reaches a magnitude such that h/q 
is no longer large compared to the nucleus, the spatial 
distribution of the nuclear spin must be considered. This 
is described by the spin structure function 



S{q) = alSooiq) + aoaiSoiiq) + alSii{q), 



(3) 



here written by using the isospin convention in terms 
of the isoscalar (oo = ap + an) and isovector (ai = 
ap — an) coupling constants. The independent form fac- 
tors, namely the pure isoscalar term Sqo, the pure isovec- 
tor term Sn and the interference term 6*01 can be ob- 
tained from detailed nuclear shell model calculations. 
Ressel and Dean [10] parameterize the structure func- 
tion in terms of y — {qh/2K)^ , where b parameterizes the 
nuclear size with b = v^fm ~2.3fm for heavy nuclei. 
For a WIMP mass of 100 GeV/c^ with velocity lO^^c, 
typical Xe nuclear recoil energies are below '--^25keV. We 
have restricted our search to below 26.9 keV, resulting 
in a maximum momentum transfer of g ~81 MeV/c, i.e 
h/q >2.4fm and y<0.25. For our analysis, we use the 
structure function calculated with the Bonn A and Ni- 
jmegen II potentials [l3| and with the QTDA method 
[13j for ^^^Xe and ^'^^Xe, respectively. 

In order to set limits on spin-dependent WIMP nu- 
cleon couplings, we follow the procedure described in 
0; 1 avoiding model-dependent assumptions on 

the (Majorana) WIMP composition. We first present ex- 
clusion limits for the cases of pure-proton (a,i ~ 0) and 
pure-neutron (op — 0) couplings, by assuming that the to- 
tal cross section is dominated by the proton and neutron 
contributions only. We calculate the 90% C.L. exclusion 
limits as a function of WIMP mass with Yellin's Max- 
imal Gap method p^ . The exclusion limits presented 
here assume a flat 19% Ceff for the xenon scintillation 
efficiency of nuclear recoils relative to electron recoils. 
The motivation for this choice is explained in reference 
[sj, in which we also show that the uncertainty in Ceff 
could raise the limits by about 15% (18%) for a WIMP 
mass of30GeV/c2 (lOOGeV/c^). 

In FiglT] we show the combined upper limit curves 
in the WIMP-nucleon cross section versus WIMP mass 
plane, for the simple case of pure neutron (left) and 
pure proton (right) couplings for the ^^^Xe and ^^^Xe 



isotopes. We also include limits from CDMS 2C|, 
ZEPLIN-II 21], KIMS [H, NAIAD [23^, PICASSO [H, 
COUPP [25] as well as the indirect detection limits from 
Super-Kamiokande [3l| for the case of pure proton cou- 
plings. Given that ^^^Xe and ^^^Xe both contain an un- 
paired neutron, XENONIO is mostly sensitive to WIMP- 
neutron spin-dependent couplings and excludes previ- 
ously unexplored regions of parameter space. The min- 



achieved at a WIMP mass of around 30GeV/c^. The 
sensitivity to pure proton-couplings is less strong, how- 
ever XENONIO improves upon the parameter space con- 
strained by the ZEPLIN-II [2l| and the CDMS ^ ex- 
periments and approaches the sensitivity of other direct 
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FIG. 1: XENONIO combined 90% CL exclusion limits for ^^^Xe and ^^^Xe for pure neutron (left) and pure proton (right) 
couplings (solid curves). The dashed curves show the combined Xe limits using the alternate form factor. Also shown are the 
results from the CDMS experiment ^ (diamonds), ZEPLIN-II [2l[ (circles), KIMS i2.2j (triangles), NAIAD 23] (squares), 
PICASSO _24j (stars), COUPP [2^ (pluses), SuperK 31] (crosses), as well as the DAMA evidence region under the assumption 
of standard WIMP nuclear recoils and dark halo parameters (green filled region) uM- The theoretical regions (blue filled) for 
the neutralino (in the constrained minimal supersymmetric model) are taken from |32|. 



detection experiments such NAIAD PICASSO 
KIMS [111 and COUPP |25|], all containing nuclei with un- 
paired protons. We also show the DAMA evidence region 
for standard WIMP interactions and halo parameters 
[m , as well as predictions of neutralino-nucleon cross sec- 
tions in the constrained minimal super-symmetric model 
(CMSSM) [s^l • Although the expected cross sections are 
still below the current experimental spin-dependent sen- 
sitivity, direct WIMP detection experiments are now for 
the first time approaching the theoretically predicted pa- 
rameter space for neutralinos. 

As a further benchmark, we consider heavy Majorana 
neutrinos with standard weak interactions. Such neutri- 
nos, with masses in the region 100-500 GeV/c^ have re- 
cently been proposed as dark matter candidates in mini- 
mal technicolor theories in cosmologies with a dynamical 
dark energy term [13, • The expected cross section on 
protons and neutrons can be written as [2^ : 



O'uN — — -J-t-l-^Cspin,!^ (4) 
TTAl 

where /i is the neutrino- nucleon reduced mass, and the 
spin enhancement factor in this case is: 

Cspin,!^ = [o-piSp) +ari('S'ri)]^ — — (5) 

with the values for the WIMP-nucleon spin factors Op 
= 0.46 and a„ = 0.34 taking into account the strange 
quark contribution to the nucleon spin, as measured by 
the EMC collaboration and given in [26| for coupling to 
protons and neutrons, respectively. 



In Fig. [2] (left) we show the predicted number of 
events in XENONIO as a function of the heavy Majo- 
rana neutrino mass, the light shaded area showing the 
excluded mass region at 90% CL for using the main 
form factors. Our result excludes a heavy Majorana neu- 
trino as a dark matter candidate with a mass between 
9.4GeV/c2-2.2TeV/c2 (9.6 GeV/c^-l.S TeV/c^ for the 
alternate form factor). We note that a heavy Majorana 
neutrino with a mass below half the Z-boson mass has 
already been excluded at LEP (30j . 

We now present the results in terms of the more gen- 
eral phase space for a„ and Op for a fixed WIMP mass. 
We follow [iSl and express the expected number of recoil 
events N^^g^o^ as a function of a„ and Opi 

Nxenon = Aa^p + Ba^an -f Ca^ (6) 

with A, B, C being constants of integration of the differ- 
ential event rate dR/dE over the relevant energy region, 
in our case 4.5keV-27keV nuclear recoil energy. 

Fig. [H (right) shows the allowed regions at 90% 
CL in the Cp — a„ parameter space for a WIMP mass 
of 50GeV/c2. We include the published CDMS [H, 
ZEPLIN-II [21I, KIMS [2^ and the DAMA allowed re- 
gion [18] for comparison. The advantage of using differ- 
ent isotopes with spin as dark matter targets is evident: 
the presence of both odd-neutron and odd-proton num- 
ber isotopes breaks the degeneracy and only the common 
space inside the ellipses is allowed by the data. 
In conclusion, we have obtained new limits on the spin- 
dependent WIMP-nucleon cross section by operating a 
liquid-gas xenon time projection chamber at LNGS, in 
WIMP search mode for 58.6 live days with a fiducial mass 
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FIG. 2: Left: Predicted number of events in XENONIO for a heavy Majorana neutrino with standard weak interaction as a 
function of the neutrino mass, using the main (sohd curve) and alternate (dashed curve) form factors. The hght shaded area 
shows the excluded mass region at 90% CL, calculated with Yellin's Maximal Gap method [l^ for the main form factors. Right: 
Regions allowed at the 90% CL in a„ — Op parameter space for a WIMP mass of 50GeV/c^ The combined limit from ^^^Xe 
and ^^^Xe is shown as a dark solid curve (using the main form factor, see text). The exteriors of the corresponding ellipses are 
excluded, the common space inside the ellipses being allowed by the data. We also show the results obtained by KIMS [23 | 
(dot-dashed), COUPP ^ (dotted) ('horizontal' ellipses) and CDMS ^ (dotted), ZEPLIN-II f?! (dashed) and the DAMA 
evidence region [l^ (light filled region) ('vertical' ellipses). 



of 5.4kg. The results for pure neutron couplings are the 
world's most stringent to date, reaching a minimum cross 
section of bxlQ-^^cm'^ at a WIMP mass of SOGeV/c^. 
We exclude new regions in the Up — an parameter space, 
and, for the first time, we directly probe a heavy Majo- 
rana neutrino as a dark matter candidate. Our obser- 
vations exclude a heavy Majorana neutrino with a mass 
between -10 GeV/c2-2 TeV/c^ for using a local WIMP 
density of 0.3GeV/cm'^ and a Maxwell-Boltzmann ve- 
locity distribution. We note that our sensitivity to axial- 
vector couplings could be strongly improved by using a 



larger mass of enriched ^^^Xe as the dark matter target. 
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